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I. INTRODUCTION 



The integrability condition called sliape-invariance originates in supersymmetric quan- 
tum mechanics The separable positive-definite Hamiltonian Hi = A is called shape- 
invariant if the condition 

A{ai)A\ai) = A^{a2)Aia2) + i?(ai) , (1.1) 

is satisfied P|. In this equation ai and 02 represent parameters of the Hamiltonian. The 
parameter 02 is a function of ai and the remainder R{ai) is independent of the dynamical 
variables such as position and momentum. Even though not all exactly-solvable problems 
are shape- invariant shape invariance, especially in its algebraic formulation has 
proven to be a powerful technique to study exactly-solvable systems. 

In a previous paper we used shape-invariance to calculate the energy eigenvalues and 
eigenf unctions for the Hamiltonian 

H = A^A + ^ [A, A^] (ag + 1) + Vm (a+A + a_A^) , (1.2) 
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where 

^± = \ (^1 ^^^2) , (1-3) 

and with i = 1, 2, and 3, are the Pauli matrices. 

This is a generalization of the Jaynes-Cummings Hamiltonian 0. A different, but re- 
lated problem was considered in Ref. |T^. Our goal in this paper is to study a further 
generalization of the Jaynes-Cummings Hamiltonian by introducing a term proportional to 
o"3 with an arbitrary coefficient (the so-called non- resonant limit). In addition to the energy 
levels we study the time-evolution and the population inversion factor. 

Introducing the similarity transformation that replaces oi with 02 in a given operator 



and the operators 



f(ai)0(ai)ft(ai) =0(02) (1.4) 



B+ = A^{ai)f{ai) (1.5) 



B_ = Bi=f\ai)A{ai), (1.6) 
the condition of Eq. ( |1 . 1| ) can be written as a commutator ^ 

B+] = rt(ai)i?(ai)f (ai) = i?(ao) , (1.7) 

where we used the identity 

Rian) = f{ai) R{an-i) f\ai) , (1.8) 

valid for any n. The ground state of the Hamiltonian Hi = A^A = B+B^ satisfies the 
condition 
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i I ^o) = = 5_ I ^o) ; 

and the unnormahzed n-th excited state is given by 
with the eigenvalue 



:i.9) 



;i.io) 



k=l 

We note that the Hamiltonian of Eq. ( |1.2| ) can also be written as 



H 



"f 


■ 




■ft 


" 


.0 


±1. 




. 


±1 



where 



h± = B+B_ + -R{ao) (ag + 1) ± VM (a+B_ + a_B. 



(1.11) 



:i.l2) 



;i.l3) 



II. THE GENERALIZED NON-RESONANT JAYNES-CUMMINGS 

HAMILTONIAN 



The standard Jaynes-Cummings model, normally used in quantum optics, idealizes the 
interaction of matter with electromagnetic radiation by a simple Hamiltonian of a two- 
level atom coupled to a single bosonic mode |1TT|-|1^. This Hamiltonian has a fundamental 
importance to the field of quantum optics and it is a central ingredient in the quantized 
description of any optical system involving the interaction between light and atoms. The 
Jaynes-Cummings Hamiltonian defines a molecule, a composite system formed from the 
coupling of a two-state system and a quantized harmonic oscillator. In this case, its non- 
resonant expression can be written as 



H = A^A + - [A, A^\ ((T3 + 1) + (<T+i + a_A^) +hAa^, 



(2.1) 



where f2 is a constant related with the coupling strength and A is a constant related with 
the detuning of the system. 

Following Ref. |p[ we introduce the operator 



(2.2) 



where the operators A and satisfy the shape invariance condition of Eq. (|1.1|) . Using this 
definition we can decompose the non-resonant Jaynes-Cummings Hamiltonian in the form 



H = Ho + H 



int 1 



(2.3) 



where 
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Ho = , 

iiint = VhQ S + ^ A (T3 . 



(2.4a) 
(2.4b) 



First, we search for the eigenstates of S^. In this case it is more convenient to work with its 
5-operator expression, which can be written as [B] 



(2.5) 



where H2 = TB-B+T"^ ■ Note the freedom of sign choice in Eq. ( p.5|) , which results in two 
possible decompositions of S^. Next, we introduce the states 



"T 


■ 




.0 


±1. 





■ 




■ft 


■ 




"^2 


" 


B+I3_ . 




. 


±1. 




. 





vl>(±)) 



T 
±1 



(2.6) 



where = C^l^[R{ai), R{a2), R{as), . . .] are auxihary coefficients and, | m) and | n) 

are the abbreviated notation for the states | ipm) and | ipn) of Eq. ( |1.1(J| ). Using Eqs. (p..7D, 
( p.5| ) and ( p.6|) , the commutation between Hi and a function of R{ak), and the T-operator 
unitary condition, one gets 



S2 I ^(±)) 



"T 


■ 




.0 


±1. 




"f 


■ 




.0 


±1. 





B+B_ + R{ao) 1 rCW | m) 
J [ C^±) I n) 

Sm + R{ao) OirCW|m) 
^nJ[cW|n) 



(2.7) 



However, using Eqs. ( |1.8| ) and ( |1.11| ) one can write 

f [S^ + R{ao)] f^ = f [R{ai) + R{a2) + ■■■ + R{aJ + R{ao)] 



R{a2) + R{az) H h R{am+i) + R{ai) = £m+i 



Hence the states 



ml 



T 
±1 



&^l\m + l)\ 



m = 0,l,2,--- 



are the normalized eigenstates of the operator 

S2 I =£ 



(2.8) 



(2.9) 



(2.10) 



We observe that the orthonormality of the wavefunctions imply in the following relations 
among the C's: 



(±) 



m+l 

/vi/(=F) I xr/(±)\ = r(±)r(=F) - r('F) = n 

\ m I m / ^m+l^m+1 " 



1 



\m Im/ ' 



(2.11a) 
(2.11b) 



Since and commute then it is possible to find a common set of eigenstates. We 
can use this fact to determine the eigenvalues of tlint and the relations among the C's 
coefficients. For that we need to calculate 



4 



H 



int 



ml \ m I 



(2.12) 



where A^-* are the eigenvalues to be determined. Using Eqs. (p.2|), ( |2.4| ) and (pl9|), the last 
eigenvalue equation can be rewritten in a matrix form as 



a 



(3 TB_ 
B+f^ -13 



T 
±1 



I m) 
|m+l) 



I m) 
|m + l) 



(2.13) 



where a = y/Ml and (3 = hA/a. Since the C's coefficients commute with the A or A' 
operators, then the last matrix equation permits to obtain the following equations 



(3 - A(^)] (f C(^±)f f\m)±a d^liTB^ | m + 1) = 
a {fC^^^f^) B+\m)T + A^f^^] Ci^i | m + 1) = . 



Introducing the operator 



B^B_ 



-1/2 



B, 



one can write the normalized eigenstate of Hi as 

I m) = (Q^Y I 0) 
and, with Eqs. ( |2.15| ) and ( p.l6| ) we can show that [ 



B+\m) = JSm.+i I m + 1) 



T5_ I m + 1) = T \m) . 

Using Eqs. ( |2.17| ), then Eqs. ( p.l4| ) take the form 



a/3 - AW] (f CWf t) ± a^E^iCi^l f | m) = 



(T'C'i^^rt) T [aP + AW] I m + 1) = 



From Eqs. (|2A8|) it follows that 



A(^) = ±aJW + /5^ 



and 



(±) 

m+l 



m+l 



(TcWf t) . 



Eqs. (imp and ( p:^ imply that 



(2.14a) 
(2.14b) 



(2.15) 

(2.16) 

(2.17a) 
(2.17b) 

(2.18a) 
(2.18b) 

(2.19) 
(2.20) 



(2.21) 
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and the eigenstates and eigenvalues of the generahzed non-resonant Jaynes-Cummings 
Hamiltonians can be written as 



and 



1 ^W) = 

1 m / 


"T 


■ 




.0 


±1. 





I m) 
|m+l) 



m = 0,1,2, ••■ 



(2.22) 
(2.23) 



a) The Resonant Limit 



From these general results we can verify two important and simple limiting cases. The 
first one corresponds to the resonant situation, for which A = (/? = 0). Using these 
conditions in Eqs. ( |2.20| ) and ( p.22| ) and Eqs. ( |2.11| ) we get 



■"m+\ 



and 



Cm+l — ^'C'm ^T''^ — Cm ~ ' 



Therefore the Jaynes-Cummings resonant eigenstate is given by 



m I 



T 
±1 



I m) 
m + 1) 



m = 0,1,2, 



(2.24) 
(2.25) 

(2.26) 



These particular results are shown in the Ref. [^. 
b) The Standard Jaynes-Cummings Limit 

The second important limit corresponds to the standard Jaynes-Cummings Hamiltonian, 
related with the harmonic oscillator system. In this limit we have that T = — > 1, 

— > a, -B+ — >■ a\ A = uj — ujo and Sm+i = {m + l)huj. Using these conditions in the 
Eqs. (ICTI ), (|]2|) and Eqs. (pi]) we conclude that 



and 



1 



a 



(±) 

m+l 



(±)^(±) ^ ^(T) 



where 



1 + 52 T 5„ 



m 

CiJ — iOn 



(2.27) 
(2.28) 

(2.29a) 
(2.29b) 



(m -|- l)r2ci; 

Therefore the standard Jaynes-Cummings eigenstate, written in a matrix form, is given by 



m I 



1 



"1 ■ 




m) 


.0 ±7^^^ 




_ m + 1) _ 



m = 0, 1,2,--- 



(2.30) 



These results are shown in many papers, in particular, in the Ref. [17 



6 



III. THE TIME EVOLUTION OF THE SYSTEM 



To study the time-dependent Schrodinger equation for a Jaynes-Cummings system in 
non-resonant situation 



d 



we can write the wavefunction as 

I = exp {~tUot/n) I 



(3.1) 



(3.2) 



and, by substituting this into Schrodinger equation and taking into account the commutation 
property between Ho and Hmf, we obtain 



ih^^ I *,(t)) = Hint I ^iit)) . 

We introduce the evolution matrix \Ji(t, 0): 

I ^,(t)) = U,(t,0) I ^,(0)). 

which satisfies the equation 



that is, in matrix form, written as 
ih 



JJ' 11' 



a 



13 TB^ 





'Un 








U22. 



(3.3) 



(3.4) 



(3.5) 



(3.6) 



where the primes denote the time derivative. One fast way to diagonahze the evolution 
matrix differential equation is by differentiating Eq. ( ^.51 ) with respect to time. We find 



in^ij.,{t, 0) = H,„i ^u,(t, 0) = ^HLu,(t, 0) 



which can be written as 



U" U" 

fj" fj" 
^21 ^22 



uji 

UJ2 



Un U12 
U21 U22 



where 



hCji = a^TB^B+T^ + (3^ = \JnVt H2 + (/lA) 



hCj2 = a^B+B_ + /32 



(3.7) 



(3.8) 



(3.9a) 
(3.9b) 



Since by initial conditions Ui(0,0) = I, then we can write the solution of the evolution 
matrix differential equation (p.7|) as 
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Ui(t,0) 



COS (cjit) sin (uit) C 
sin (cj2t) D cos {uJ2t) 



and the (7 and D operators can be determined by the unitarity conditions 

u|(t,o)u,(t,o) = u,(t,o)ul(t,o) = i. 



In the appendix A we show that the unitarity conditions ( 3.11 ) imply 

i 



C = -D^ 



{H2 



1/4 



(3.10) 

(3.11) 

(3.12a) 
(3.12b) 



D = -CK 

Therefore, we can write the final expression of the time evolution matrix of the system as 



U,(t,0) 



cos {Coit) sin {uJit) C 
— sin {uj2t) cos (1^2^) 



(3.13) 



For Jaynes-Cummings systems an important physical quantity to see how the system 
under consideration evolves in time is the population inversion factor [|ri|,|l^,|l^ , defined by 

W(t) = a+{t) a^it) - a4t) a+(t) = asit) , (3.14) 

where the time dependence of the operators is related with the Heisenberg picture. In this 
case, the time evolution of the population inversion factor will be given by 



da-t(t) 1 - + , , r " 1 " / 
— ^ = -Ul t,0 a3,H U, t,0 
at in L J 



and since we have 

(T3, H = a 
then Eq. ( p.l5| ) can be written as 

dt h 



o"3, S = —2a S o"3 , 



(3.15) 



(3.16) 



(3.17) 



We can obtain a differential equation with constant coefficients for cr^it) by taking the time 
derivative of Eq. (p.l7|) 



d^a3{t) 2ia j dS{t) 



dt^ 



h 



dt 



da3{t) 
dt 



(3.18) 



Having in mind that 



ffl4uJ(,0)[S,H]u.(,0) 



(3.19) 



and that 
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we conclude that 



S,H 



dS{t) 



a[3 S,(3"3 =2q;/?S(3"3 



2ia(3 



dt h 

Using Eqs. ( |XT7| ) and into Eq. {^J^ we obtain 



where 



02 



4«2 



+ e'a3(t) =F(t) 



. 4«'/^ TTt 



uKt,o)su,(t,o). 



(3.20) 
(3.21) 

(3.22) 

(3.23a) 
(3.23b) 



Eq. ( p.22|) corresponds to a non- homogeneous hnear differential equation for a^it) with 
constant coefficients since and H commute and, therefore, is a constant of the motion. 
The general solution of this differential equation can be written as 

a3(t)=a^(t) + a^(t), (3.24) 

and each matrix element of the homogeneous solution, satisfies the differential equation 



+ z>|ag(t) = 0, j,k = l,OT2 



dt^ 



with 



hv2 = 2ayB+B^ = 2\/hQHi . 
The solution of Eq. ( p.25| ) is given by 

^fkit) = Vjit) Cjk + Zj{t) djk ■ 

where 



yj{t) = cos (uji) 
Zj(t) = sin (z>jt) , 



(3.25) 

(3.26a) 
(3.26b) 

(3.27) 



(3.28a) 
(3.28b) 



and the coefficients Cjk and djk can be determined by the initial conditions. 

The matrix elements of the particular solution of the (T3(t) differential equation need to 
satisfy 



dt^ 



+ 9^ af,,{t) = Fjk{t) , J, fc = l,or2 
9 



(3.29) 



and can be obtained by the variation of parameter or by Green function methods, giving 

^fkit) = ^7' {%W l^d^m)FAO-ml^d^^,iOFjkiO} , (3.30) 

where we used that the Wronskian of the system of solutions yj{t) and Zj{t) is given by Oj. 

After we determine the elements of the F(t)-matrix, it is necessary to resolve the integrals 
in Eq. ( ^.301 ) to obtain the explicit expression of the particular solution. In the appendix 
B we show that, using Eqs. (p.2|), (|3.13|) , and ( p.23|) , it is possible to conclude that these 
matrix elements can be written as 



2 



+ ^1^1 ^hI''^ [^ziit) Qscit' '>i,u}i,u}2) - yi{t) Qccit'^ 01,^1,^2)} \l B+f^ , (3.31a) 
^uit) = ^-K^\lf¥^{z2{t) Qccit-^ y2.Cj2.Cjx) - Ut) g^scit; i>2,Cj2,Cj,)} ^/fH 

+ |z>r^i^2^' {zi(t) gi'Jit; h,C;,,u;2) + yi{t) g'bkt-^ vi,Cji,Cj2)]hI'^ , (3.31b) 



^2l{t) = |z>2-V^+^^ {Zl{t) g^^^t- KCJi,CJ2) - yi{t) gicit; h,Cjl,U2)} V^+f t 

+ ^K^hI^^ {z2{t) gssit; j)2,Cj2,Cji) - y2{t) ^{tsH^; i'2,uj2,Cji)} H^^^ , (3.31c) 

'^2''2W = ^1^2"'/^ {^iW gcsit; h,C;i,C;2) - yi{t) g^ss\t; h, Cj^ UJ2)} hI^' 

+ i^V2^Hl'^ [hit) gi'ckt; h,C^2,C^i) + Ut) gccit; U2,Cj2,Cj,)} ^/fE , (3.31d) 

where 7 = Aa^jS/fi^^ and the auxiliary functions are given by 

g^}kt-p,q,r)=TxY{t\p-q,r)±TxY{t\P + q,r), X,Y = C 01 S , (3.32) 

with 

J-'ccit] x,w) = cos(x^) cos(w^) 

Jo 



^2m^2n ^2m+2n+l 



Y.J l)"^^"(2,„),(2n)!(2m + 2n + l) ^^'^^^^ 

ft 



J-'csit; x,w) = cos(x^) sin(w^) 

Jo 



00 ^2m,",2n+l j-2m+2n+2 

y (-1)"^+" ^^^^ --^ (3.33b) 

„tlo (2m)!(2n+l)!(2m + 2n + 2) ^ ^ 



^sc{t] x,w) = f d^ sin [x^) cos (w^) 
Jo 
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00 -2m+l o 2n i2m+2n+2 

m^J ' (2m+l)!(2n)!(2m + 2n + 2) > 

rt 



^ss{t]x,w) = [ d^ sin(a;^) sin (w^) 
Jo 



00 -2m+l o 2n+l j.2m+2n+3 

y J!i ^ ^. (3.33d) 

mn=o (2m+l)!(2n + l)!(2m + 2n + 3) ^ ^ 
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With these results for the particular solution we can conclude that 

aj(0) = 



dt 



Now, using Eqs. (|3.17|) , (|3.24|) , (|3.27|) , (|3.34|) and the initial conditions, we have 



'^3(0)] 



ij ' '^ij 



dt 



2ia 
IT 



S(0)a3(0) 



z/j dij . 



(3.34) 

(3.35a) 
(3.35b) 



Therefore, the final expression for the elements of the population inversion matrix of the 
system can be written as 



2ia 



[a^m, = cos (z>,t) [^3(0)]., + — sin (z>,t) z/r^ S(0) ^3(0) . . + afJt) . 



(3.36) 



Again, using these final results we can verify two important and simple limiting cases, 
a) The Resonant Limit 

The first one corresponds to the resonant situation (A = 0). Eqs. ( |3l9|) , (|3.13|) , (|3.26|) 
and ( p.31| ) allow us to conclude that, in this case, the evolution matrix of the system is given 
by 



Ui(t,0) 



cos 
sin 



(|z>it) sin(iz>it) 



(iz>2t) cos (|z>2t 

and the elements of the population inversion of the system are 

2ia r- 
[a3{t)]ij = cos (Oit) MO)],^ + — sin (z>,t) Or^ [S(0) ^3(0) 



(3.37) 



(3.38) 



b) The Standard Jaynes-Cummings Limit 

This second important limit corresponds to the case of the harmonic oscillator system, 
and in this limit we have that T = T'^ — > 1, — > d, -B+ — > a) and [a, d^] = huj. 
With these conditions the operators uoi and Cj2 commute, and this fact permits to evaluate 
the integrals related with the particular solution of the population inversion elements using 
trigonometric product relations. Using that and the expressions obtained in the appendix 
B, after a considerable amount of algebra and trigonometric product relations we can show 
that is possible to write the expressions for the (5"j^(t)-matrix elements as 



^fi(^) = ^2^1 ^Vd{ICs{t;u2,i^i,i>2) - }Csit;u2, -^1,1)2)} (da^) 

- i^i/f ^ (dd^) ^ {ICs{t-,uj2,uJi,i'i) - l<^s{t]'^2,~i^i,i>i)}^ (3.39a) 
^uit) = ^i>i^Vd{ICcit;uj2,u}i,i'2) - }Cc{t; LJ2, -toi, h)} Vd 
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^z>i ^ (aa"^) {/Cc(t; t^i, z^i) - /Cc(t; uj2, -t^i, ^^i)} (a'^a) 



— — i>n ""^ I a'a 



1 /;it; 



1/4 



{fCcit; Qj2, uji, z>2) - /Cc(t; c2'2, -t^i, '>2)} (aa"^) 



1/4 



1/4 



^22(^) = K^^2 {^sit] ^2, t^i, z>i) + KLsit] ^2, -t^i, i^i)} [a^aj 



1/4 



- i^z>2 ^ {a)aY'^ {lCs{t\uJ2,(^i,V2) + K,s{t\ uj2, -uji, 1)2)} , 

where, now, the auxihary functions are given by 

f sin [{j) + q)t\ — {p + q) sin {ft) 



ICs{t;p,q,f) 



)Cc{t;p,q,f) 



f2 — (]3 + q) 
r cos [{p + q)t] — f cos (ft) 
f2 — (p + g)^ 



(3.39b) 
(3.39c) 
(3.39d) 

(3.40a) 
(3.40b) 



Considering the expressions above we may easily verify that the particular solution for the 
population inversion factor must still satisfy the initial conditions ( p.34| ). Therefore, in this 
case the final expression for the population inversion factor has the same form given by 
Eq. i ^Ml . with 



tiC'2 = 2aV d^a . 



huoi = a J da) + 



huj2 = OL\l d)d + . 



(3.41a) 
(3.41b) 



IV. CONCLUSIONS 

In this article we extended our earlier work on bound-state problems which rep- 
resent two-level systems. The corresponding coupled-channel Hamiltonians generalize the 
Jaynes-Cummings non-resonant Hamiltonian. If we take the starting Hamiltonian to be the 
simplest shape-invariant system, namely the harmonic oscillator, our results reduce to those 
of the standard non-resonant Jaynes-Cummings approach, which has been extensively used 
to model a two-level atom-single field mode interaction whose detuning it is not null. 

Another possible extension of our model is to consider intensity-dependent couplings. 



This will be taken up in the following paper |T8 
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Appendix A 

Here we give the steps used to obtain the specific form of the operators C and D. Using 
Eq. (|3.10|) into the unitary condition equation actually we can show that the C and 



D operators need to satisfy the following six conditions 



CC^ = C^C = 1 (la) 

DD' = D^D = 1 (lb) 

D^^ sin {LJ2t) = — sin (uit) C (Ic) 

b cos {Cjit) = - cos {uj2t) . (Id) 



At this point we can use the following property 



TB_ U2 = \JTB_ ^a'^B+B_ + 

a^fB_B+B_ + fB_(3^/h 



a^TB.B+T^TB^ + p'^TB./h 



a^TB^B+T^ + P'^/h \JTB_ 



ooi ^TB_ . (2) 



Then, with this result we have 



ITB_ CjI = ^TB_ U02 Qj2 = u^i yTB_ UO2 = u)\ ^TB_ , (3) 
and finally, by induction, we conclude that 



lTB_ 00^ = cOi" \JTB^ . (4) 
In the same way. 



uji = \JB+T^ ya^TB^B+T^ + 

a^B+f^fB_B+f^ + B+f^pyn 
a^B+B_B+f^ + (3^B+fyn 



0j2 \JB+T\ . (5) 
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Then, with this result we have 



'B+T^ CjI = \/B+T^ LJi ui = UJ2 \/B+T^ LJi = CjI ^B+T^ , (6) 
and finally, again by induction, we get 



fi+Tt Co'l = 00^ \JB+T^ . (7) 

Using the properties given by Eqs. (H) and (|^) and the forms of C, D operators, defined by 
Eqs. ( |3.12| ), we can verify that 

-i) 1 /TT 1 



CC^ = D^D = -^JTB^ ^B^T^^^ = ^V^2^ = 1 , (8) 

and 



^1/4 V - V + ^1/4 ^1/4 V .^1/4 



C^C = DD^ = V^yi iui^- = ^2^^ = 1 . (9) 

Also using the series expansion of the trigonometric functions, we can show that 

[Uj2t) 



i,tsi„(^.o-^v'™-i:(-i)'p„^„. 



00 



UJit) —I 



= - sin {uJit) C , (10) 

where we used the commutation between and Cj\ (see Appendix B). In the same way 
we can prove that 



b cos(c^,t) = v^7^E(-ir^ 



2n 



Y.{-lr^JB^T^ 



(2n)! Hll^ 



n=0 v-^'t.;. 
- COS {ib2t) . (11) 



Again, we used the commutation between H2 and a)i. 
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Appendix B 



In this appendix we sliow the necessary steps to obtain the explicit expressions of the 
particular solution elements of the population inversion factor. To resolve the integrals in 
Eq. (|3.30| ), first we need to determine the elements of the F(t)-matrix. To do that we can 
use Eqs. ( ^231) , and ( ^TsD to write down 
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7 |cos (uit) TB^ sin (0)2^) C*^ + (7 sin (0)2^) -B+T^ cos (ooit)^ 

TB^ cos {Cj2t) sin {ujit) H2^'^ — H2^'^ sin (uoit) cos {uj2t) \J B+T^ 



Fuit) = 7 |cos (uJit) TB^ cos {uj2t) — C sin (ci)2t) B^T"^ sin {Cjit) 

ItB^ cos {uj2t) COS {Cjit) \JtB^ + H2^ sin sin {Cj2t) H\^'^ 



7 



F2i(t) = 7 {cos (t2'2t) 5+f ^ COS pit) - (7"^ sin {(bit) fB_ sin (^2^2^) C*"^} 



7 



B^T'i cos (cuit) cos (ct;2t) \/ -B+Tt + if^'^^'^ sin (c2'2t) sin (cDit) 



ri/4 



F22(t) = 7 {(7"^ sin (uit) TB^ cos (tl^st) + cos {u2t) B+f^ sin (uit) C} 



Z7 



-B+Tt cos (ooit) sin (ci)2^) -f^i^"^ — H^'^ sin (ci)2^) cos \/TB 



1/4 



(12a) 



(12b) 



(12c) 



(12d) 



where 7 = Aa^jS/h^. Here we used the properties (|I]), (^ and (|^, together with the following 
operators relations 



CxB^T^ 



B^TW 



fB_ = iHI'^ 



(13a) 
(13b) 



Now, keeping in mind that cuj] = 0, (j = 1, or 2), so we may use the trigonometric rela- 
tionships involving products of trigonometric functions with arguments Vjt and ujjt (since we 
have exp {ujt) exp (±ci>jt) = exp [{uj ± Cjj)t\). Then, using those relationships, the following 
commutators 



V2,Hi = U2,H 







(14) 



and the same properties (|l]), (H) and (0), we can show that 

yiit)Fii{t) 



.7 



i-\jTB^{co^[{v2- Cj2)t\ sin (cJit) + COS [(z>2 + cj2)t] sin (c2;it) } if2 



1/4 



+ i^H^^ {sin [{vi — Cji)t] cos {tli2t) — sin [(z>i + Coijt] cos (0)2^)} \l -B+Tt (15a) 

yi{t) Fi2{t) = ^\JfB_ {cos [(z/2 - ^^2)^] COS {Cjit) + cos [(Z>2 + CJ2)t] COS (a)it)} \JfB_ 

+ ^^2'''^ {sin [('^i + (Ji)t\ sin (ct;2t) — sin [(z>i — cji)i(:] sin (1^2^)} H\^'^ (15b) 

hit) F2l{t) = '^\/B+f^ {cos [(z/i - 0di)t] cos (cJ2t) + COS [(z/i + CJi)t] COS {i02t)} \J Bj^f^ 

+ ^-ff^^"^ {sin [(z>2 + c2;2)t] sin [Cjii) - sin [(z>2 - Cj2)t\ sin (cJit)} i^a^"^ (15c) 
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.7 



hit) F22{t) = i-y B+T^ {cos [(z>i - Qji)t\ sin {Cj2t) + cos [{vi + Cji)t\ sin (^^2^)} -^i 



1/4 



+ z-^i^^y^ {sin [(z/2 — Cj2)'t] cos (ci)it) — sin [(i>2 + '^'2)^] cos (cuit)} \/ Ti?_ . (15d) 



In a similar way, we can show that 



.7 



-Fii(t) = i-\jTB_ {sin [(z>2 - cj2)t] sin (c^it) + sin [(z>2 + cj2)t] sin [Coit)] H2 



1/4 



— z-^i^2^^ {cos [(z>i — ci)i)t] COS (0)2^) — COS [(z>i + iUi)t] COS (ii'2i)} \/ B^T'f (16a) 
A2(t) = {sin [(z/2 - cl's)^] cos (cl^it) + sin [(z>2 + ^^2)^] cos {Cjit)}^fl3^ 



l^H^'^ {cos [(z>i + sin (cj2^) — cos [(z>i — c2'i)t] sin {Cj2t)} H\^'^ 

7 



(16b) 



^2(t) F2i{t) = ^^B+T^ {sin [(z>i - ui)t] cos (cl;2t) + sin [{Oi + ui)t] cos (^2^2^)} ^ B+T^ 

— {cos [(z/2 + ^^2)*^] sin (cjit) — cos [(z/2 — a>2)t] sin [Cjit)} Hy'^ (16c) 



52(t) -^22(i) = i^\J B+f'^ {sin [(i>i - uJi)t] sin (cj2t) + sin [(i>i + u)i)t] sin (0)2^)} -f^i 



1/4 



— i^H\^^ {cos [(z/2 — ci)2)t] COS iCjit) — COS [(z/2 + 0)2)^] COS (ci)it)} \JtB_ . (16d) 

The non-commutativity between the operators uJi and CD2 imply that to calculate the 
integrals involving the terms given by Eqs. (|TB|) and (|TB|) we need to use the series expansion 
of the trigonometric functions. In this case the integrals can be easily done because the 
time variable can be considered as a parameter factor. Finally, inserting these results into 
Eq. ( p.30|) it is trivial to find the expression ( p.31|) for the matrix elements of the particular 
solution. 
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